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ABSTRAGT 

We give an overview of the Grism Lens Amplified Survey from Space (GLASS), a large Hubble 
Space Telescope program aimed at obtaining grism spectroscopy of the fields of ten massive clusters 
of galaxies at redshift z = 0.308 — 0.686, including the Hubble Frontier Fields (HFF). The Wide Field 
Gamera 3 yields near infrared spectra of the cluster cores, covering the wavelength range 0.81 — I.69//m 
through grisms G102 and GI4I, while the Advanced Gamera for Surveys in parallel mode provides 
G800L spectra of the infall regions of the clusters. The WFG3 spectra are taken at two almost 
orthogonal position angles in order to minimize the effects of confusion. After summarizing the 
scientific drivers of GLASS, we describe the sample selection as well as the observing strategy and data 
processing pipeline. We then utilize MAGSJ07I7.5+3745, a HFF cluster and the first one observed by 
GLASS, to illustrate the data quality and the high-level data products. Each spectrum brighter than 
Hab = 23 is visually inspected by at least two co-authors and a redshift is measured when sufficient 
information is present in the spectra. Furthermore, we conducted a thorough search for emission lines 
through all the GLASS WFG3 spectra with the aim of measuring redshifts for sources with continuum 
fainter than Hab = 23. We provide a catalog of 139 emission-line based spectroscopic redshifts for 
extragalactic sources, including three new redshifts of multiple image systems (one probable, two 
tentative). In addition to the data itself we also release software tools that are helpful to navigate the 
data. 

Subject headings: gravitational lensing: strong 


I. INTRODUCTION 

Over the past 25 years, the Hubble Space Telescope 
(HST) has assembled a phenomenal legacy of extragalac¬ 
tic surveys. Several legacy fields have been imaged at 
a variety of wavelengths, becoming the focus of ground 
based telescopes and communities interested in interme¬ 
diate and high redshift science (Ferguson et al. 2000; Gi- 
avalisco et al. 2004; Koekemoer et al. 2011; Scoville et al. 
9; Ellis et al. 2013). Additionally, pure-parallel fields 
surveys have covered large uncorrelated areas providing 
even larger samples and means to control sample vari¬ 
ance (Trenti et al. 2011; Schmidt et al. 2014a; Golbert 
et al. 2013). 

A perhaps less well-known part of the HST legacy are 
spectroscopic surveys carried out initially with the Near 
Infrared Camera and Multi-Object Spectrograph (NIC- 
MOS McCarthy et al. 1999) and the Advanced Camera 
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for Suverys (ACS Pirzkal et al. 2013) and more recently 
with the Wide Field Camera 3 (WFC3 Atek et al. 2010; 
Brammer et al. 2012; Morris et al. 2015). By virtue of 
its low background and exquisite image quality compared 
to the ground, HST excels in low resolution slitless spa¬ 
tially resolved spectroscopy, particularly in the near in¬ 
frared. Thus, it has become a workhorse for galaxy evo¬ 
lution studies, including: the study of spatially resolved 
star formation (Nelson et al. 2012); the detection of rare 
galaxies with strong emission lines and faint continuum 
emission (Atek et al. 2011); and the measurement of ab¬ 
sorption line redshifts for z > 1 galaxies (van Dokkum 
et al. 2011; Newman et al. 2013) 

Current HST spectroscopic surveys have focused on 
relatively shallow spectroscopy of extragalactic fields. 
Although this is a very productive way to survey large 
areas of the sky, the target depth has typically been insuf¬ 
ficient to address some of the most interesting questions 
in galaxy formation and evolution. For example, HST 
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spectroscopic surveys have been able to detect Lya for 
the brightest galaxies usually only during the tail-end of 
cosmic reionization (Rhoads et al. 2013). Additionally, 
at lower redshifts, legacy fields and random parallels do 
not contain massive clusters of galaxies and thus cannot 
probe the densest environments. 

The fields of massive clusters of galaxies provide a very 
powerful complement to the traditional legacy fields. On 
the one hand they contain massive clusters of galax¬ 
ies and thus enable the study of environmental effects. 
On the other hand the clusters themselves act as grav¬ 
itational telescopes by magnifying objects in the back¬ 
ground through the gravitational leasing effect. The 
leasing magnification enables the study of intrinsically 
fainter and more compact galaxies in the background 
than would be possible in empty fields, at the expense of 
increased contamination by the foreground and a loss of 
cosmic volume for a fixed solid angle. It is important to 
stress that the net effect of a gravitational telescope on 
the number counts of a background population (some¬ 
times known as magnification bias) can be both positive 
or negative, depending on the slope of the luminosity 
function: if the differential luminosity function is steeper 
than dn/dL oc 1/“^, one detects more objects behind 
a gravitational telescope, if it shallower the empty field 
has a higher yield (e.g., Meylan et al. 2006). In prac¬ 
tice - at the bright end of the luminosity function where 
counts are rising exponentially - the net gain from lens- 
ing is huge and aids greatly in the detection of galaxies at 
z > 6. Even more importantly, for a fixed exposure time, 
lensing allows one to probe an intrinsically fainter popu¬ 
lation, thus providing qualitatively different insights into 
the source population. 

Exploiting lensing magnification (Treu 2010; Coe et al. 
2014) is one of the motivations behind many HST cluster 
surveys (e.g.. Smith et al. 2001; Postman et al. 2012; 
Bayliss et al. 2014) and the main driver of the ongoing 
Hubble Erontier Eield initiative (Lotz et al. 2015, in 
preparation). 

The Grism Lens Amplified Survey from Space 
(GLASS; PI Treu; GO 13459) has been designed and 
carried out in order to address some of the most com¬ 
pelling scientific questions by combining the power of 
deep HST spectroscopy with the magnification effect of 
gravitational lensing. The key science questions to be 
addressed by GLASS are 

1. How and when did galaxies reionize the Universe 
(if they did)? 

2. How do gas and metals cycle in and out of galaxies? 

3. How does galaxy evolution depend on characteris¬ 
tics of the local environment? 

Eirst results on these topics are given in the papers by 
Schmidt et al. (2015), Jones et al. (2015), and Vulcani et 
al. (2015, in prep). 

In addition to the key questions driving the survey de¬ 
sign, the data collected by GLASS enable the study of 
luminous and dark matter in the cluster themselves, by 
providing spectroscopic redshifts of background sources 
to improve the lensing models, and provide additional 
epochs to search for high redshift supernovae gravita¬ 
tional lensed by the clusters. Examples of these two ad¬ 
ditional science goals are given in the early science papers 



Fig. 1.— Predicted number of Lya emitters per WFC3 field in 
cluster and blank fields based on the conditional probability func¬ 
tion of Lya emission for LBGs, and the LEG luminosity function 
model by Mason et al. (2015). Measuring the evolution of the lu¬ 
minosity function of Lya emitters is one of the goals of GLASS. At 
the depth achievable with the WFG3 G102 grism lensing magnifica¬ 
tion has a positive effect at all relevant redshifts. The dashed lines 
are for blank fields, the solid lines are for specific GLASS clusters 
(A2744=black; A370=red; MAGS0416=green; MAGS0717=blue; 
MAGS1149=cyan; RXJ1347=yellow; RXJ2248=mapnta). At ^ - 
6 the conditional probability distribution function is as described 
by Treu et al. (2012), based on the data by Stark et al. (2011). At 
z 7 we use the two “patchy” measurements by Treu et al. (2012, 
upper curve of each pair) and Pentericci et al. (2014, lower curve) 
to give a sense of the measurement uncertainties. At 2 : ~ 8 the 
conditional probability has been assumed to be the same as 2 ; ~ 7. 
More details on the models are given in the main text. 

by Schmidt et al. (2014b, hereafter paper 0), Kelly et al. 
(2015), and Wang et al. (2015). 

In this paper we give an overview of GLASS 
and we present the first release of the data for 
MAGSJ0717.5+3745, the first cluster targeted by the 
survey. The paper is organized as follows. In Section 2 
we review the scientific drivers of the survey. In Section 3 
we give the rationale for selecting the sample of clusters 
and summarize their main properties. In Section 4 we 
detail the observational strategy as well as the data re¬ 
duction. In Section 5 we describe the postprocessing of 
the spectra, including procedure and tools for visual in¬ 
spection and redshift determination. Section 7 concludes 
with a summary. Present and future data products are 
and will be made public through the HST archive. All 
magnitudes are given in the AB system and a standard 
concordance cosmology with = 0.3, ft a = 0.7 and 
h = 0.7 is adopted when necessary. 

2. SGIENTIFIG DRIVERS 

2.1. Key Science Driver 1: Gas and galaxies at the 
epoch of reionization 

The Universe is known to undergo dramatic evolution 
in the redshift range ^ = 11 — 6 while becoming fully 
ionized (Bennett et al. 2013; Planck Collaboration 2015; 
Robertson et al. 2015; Mason et al. 2015). However, 
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how and when exactly it gets there is the topic of much 
debate. The interaction between the first galaxies and 
the intergalactic medium (IGM) is poorly understood, 
leaving many questions unanswered. How many ionizing 
photons escape from the local interstellar medium (ISM) 
into the IGM? What is the topology of the ISM/IGM -g 
at these redshifts? What is the relative contribution of 
galaxies vs active nuclei to ionizing budget (Madan et al. b 
2004; Ricotti & Ostriker 2004; Giallongo et al. 2015)? co 

The reduced observed Lya flux from galaxies at 2 : > 6 ^ 

indicate that the process is still ongoing dit z ^ 6 — 7 o 

(e.g., McQuinn et al. 2007; Dijkstra et al. 2011; Jensen <0 

et al. 2013; Choudhury et al. 2014). This indication for S 

“late” reionization is consistent with constraints derived 5 

from the Lya forest (e.g., Becker et al. 2015) and also z 

observed non-detection of the kSZ effect on small scales 
in the cosmic microwave background (Mesinger et al. 

2012). In detail, quantitative constraints inferred from 
Lya emitting galaxies still suffer significantly from ob¬ 
servational uncertainties (e.g., Mesinger et al. 2015). In 
addition, uncertainties exist on the theoretical side which 
include how Lya photons escape from galaxies, how this 
depends on galaxy properties, the role of galactic winds, 
and how much Lya is reprocessed by the GGM (giving 
rise to ‘blobs/halos’ see, e.g., Dijkstra 2014, for a review). 

Recent ground-based observations suggest dramatic 
and rapid evolution of the IGM/ISM beyond z ^ 6. 

An important clue is that the fraction of dropouts that 
are Lya emitters, increases steadily out to z ^ 6 (e.g.. 

Stark et al. 2011), while it declines significantly at z ^ 7 
and above (Kashikawa et al. 2006; Fontana et al. 2010; 
Pentericci et al. 2011; Schenker et al. 2012; Ono et al. 

2012; Glement et al. 2012; Treu et al. 2012, 2013). Large 
samples are now available at redshifts up to 7 where 
high-sensitivity multiplexed optical spectrographs can be 
used, and suggest that the IGM might indeed be partially 
neutral with a patchy distribution of Lya optical depth 
(Pentericci et al. 2014, hereafter P14). 

In contrast, spectroscopic follow-up at 2 ; > 7 has been 
very difficult, owing to the atmosphere which limits stud¬ 
ies to a few specific redshift windows, inspite of the recent 
introduction of high sensitivity high multiplexed near 
infrared spectrographs (e.g., MOSFIRE McLean et al. 

2012). This difficulty is compounded by the weak aver¬ 
age Lya emission for the brighter continuum sources typ¬ 
ically accessible with spectroscopy in legacy fields com¬ 
pared to the fainter ones (Stark et al. 2011). 

The goal of GLASS is to target a large sample of in¬ 
trinsically faint Lyman break galaxies (LBGs) at 2 ; > 6.5 
in order to make the most accurate measurement to date 
of the distribution of the Lya optical depth, and thus 
probe the properties of the IGM/ISM. The first result 
of this investigation are discussed in detail by Schmidt 
et al. (2015). To achieve this goal GLASS exploits leasing 
magnification by the foreground clusters to reach fainter 
luminosities than would be possible in blank fields and 
thus in principle larger number of detections. Further¬ 
more, theoretical (Zheng et al. 2011) and observational 
(Ouchi et al. 2009; Guaita et al. 2015) arguments indi¬ 
cate that high- 2 ; Lya emitters might be significantly more 
spatially extended than the associated UV continuum. 
GLASS should provide the first constraints on the spa¬ 
tial extent of the Lya emission on subarcsecond scales, 
via an individual study of magnified bright sources and 
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Fig. 2.— Number of Lya lines with flux above 3-cr as a function 
of total G102 orbits divided over clusters. Once a sufRcent 
depth is reached, the number of detection increases linearly with 
the number of clusters and less rapidly with exposure time owing to 
the loss of efficiency in magnification bias in addition to the normal 
scaling of noise and number counts observed in empty fields. The 
predictions are based on the model shown in Figure 1 and described 
in the main text. Only the results for cluster RXJ1347 are shown 
for clarity. At 2 ; ~ 7 and 8 the upper and lower lines of each color 
represent Cp = 0.66 and Cp = 0.46 respectively, as in Figure 1. 

a stacking analysis of the sample Schmidt et al. (2015) 

2 . 1 . 1 . Predicting the Lya number counts 

As a reference. Fig. 1 shows the predicted number of 
dropouts as a function of Lya flux. The forecast is based 
on three main ingredients: i) the luminosity function of 
LBGs; ii) the conditional probability of the Lya equiva¬ 
lent width for LBGs of a given continuum magnitude; iii) 
the distribution of magnifications provided by the fore¬ 
ground lensing clusters. For the first ingredient we use 
the model of the LBG LF recently proposed by Mason 
et al. (2015) building on earlier work by Trenti et al. 
(2010) and Tacchella et al. (2013). This model is in ex¬ 
cellent agreement with with the actual number of LBGs 
detected in blank fields (e.g., Bradley et al. 2012; Oesch 
et al. 2013; Schmidt et al. 2014a; Bouwens et al. 2014). 
For the second ingredient we use the parametrization of 
the distribution of Lya equivalent width introduced by 
Treu et al. (2012, hereafter T12). Based on data from 
the literature (Stark et al. 2011) at 2 : ^ 6 the intrinsic 
rest-frame distribution is: 

Pe{W) = -^^e-i(^cTH{W) + (1 - A)5{W), (1) 
y2TiWc 

with Wc= 47A, A=0.38 for sources with —21.75 < 
Muv < —20.25 and Wc=47A, A=0.89 for sources with 
—20.25 < Muv < —18.75. A is the fraction of emitters, 

S is the delta function, and H is the step function. As 
an illustration and to avoid overcluttering the figure, we 
consider here the “patchy” family of models introduce 
by T12, which appears to provide the best description of 
the data at z ^ 7. For this family of models, the intrinsic 
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distribution of equivalenth widths is expressed in terms 
of pq and the parameter ep: 

Pv{W) = e^p^{W)^{l-e^)5{W) = 

As an illustration of observational uncertainties we con¬ 
sider two values of ep = 0.66,0.46, taken from T12 and 
P14, respectively. At z ^ 8, only upper limits are avail¬ 
able (Treu et al. 2013). For simplicity we assume no 
evolution between z ^ 1 and z ^ 8, even though the 
number counts might be significantly lower at z ^ 8. In 
any case, the numbers are sufficiently large that a smaller 
number of detections at z ^ 8 would pose a stringent 
limit on the Lya optical depth. For the third ingredi¬ 
ent we use magnification maps based on grid-based weak 
and strong lensing models constructed by members of 
our team (e.g., Bradac et al. 2008; Wang et al. 2015, 
and Hoag et al. 2015, in preparation), using the code 
developed by Bradac et al. (2005), in the adaptive ver¬ 
sion (Bradac et al. 2009). Those are available for 7/10 
clusters at the moment of this writing, including both 
HFF and non-HFF clusters. These 7 clusters are rep¬ 
resentative of the full GLASS sample and therefore the 
range of magnification maps gives a good sense of range 
of expected variations from cluster to cluster. Mock cat¬ 
alogs of high-redshift galaxies are ray traced through the 
magnification maps to generate suitable mock catalogs 
in the image plane taking into account the boost in flux 
as well as the change in solid angle. The Lya emission 
is assumed to be unresolved for this excercise, consistent 
with the observations (Schmidt et al. 2015). 

The figure shows that magnification bias is positive 
in all cases, and cluster fields always yield more sources 
than blank fields, in the range of fluxes and redshifts 
considered here. Furthermore, there are cluster to clus¬ 
ter variations, but they are not extreme. This is due 
to the fact that the magnification effect of the most ex¬ 
treme lensing clusters, such as those in the HFF sample, 
is somewhat mitigated by the size of the field of view of 
WFC3. When the critical lines are larger than the WFC3 
field of view, the high magnification areas can be lost to 
a single pointing campaign. As one can see from the 
range in the predictions, cluster to cluster variations are 
comparable to current uncertainties on the distribution 
of equivalent width of Lya. 

2.2. Key Science Driver 2: The Gas and Metal Cycles 
of Galaxies 

How gas flows in and out of galaxies is a central ques¬ 
tion in galaxy formation and evolution as emphasized 
for example by the most recent Decadal Survey (Council 
2010). The GLASS survey is particularly well suited for 
addressing gas flows at redshifts z 2, corresponding 
to the peak period of cosmic star formation rate den¬ 
sity and hence rapid gas inflow and outflow rates (e.g., 
Madau & Dickinson 2014; Shapley 2011). We have iden¬ 
tified a series of specific questions that can be addressed 
by GLASS in order to better understand the role of gas 
flows. 

How do metallicity gradients evolve with cosmic time? 
Metallicity gradients are sensitive to the history of bary- 
onic assembly: gas accretion, mergers, star formation. 


and outflows. Gradients measured from spatially re¬ 
solved emission line diagnostics have provided significant 
insight into the assembly process at z = 0 (e.g., Vila- 
Costas & Edmunds 1992; Rupke et al. 2010; Bresolin 
et al. 2012). Cosmological simulations coupled with 
chemical evolution models can reproduce these data, but 
various models predict different behavior at earlier times 
(e.g., Pilkington et al. 2012; Gibson et al. 2013). Mea¬ 
surements for a large sample of galaxies over a range 
of redshift thus provide important constraints on galaxy 
assembly history (e.g., Troncoso et al. 2014). 

How does the mass-metallicity-SFR relation evolve at 
low stellar mass? Local galaxies he on a tight rela¬ 
tion where metallicity increases with stellar mass (e.g., 
Tremonti et al. 2004) and typically decreases with SFR 
(e.g., Mannucci et al. 2010). Both effects can be ex¬ 
plained by the cycling of gas and metals between galax¬ 
ies and the IGM and are more pronounced at low mass. 
The cosmic evolution of this relation at low stellar mass 
is thus a powerful test of theories of feedback and of 
metal enrichment of the IGM (e.g., Finlator & Dave 2008; 
Henry et al. 2013). 

What causes the offset of star forming galaxies in the 
BPT diagram (Baldwin et al. 1981) at high redshifts? 
Star forming galaxies and AGN at z = 0 he on separate 
well-defined loci in diagnostic line-ratio diagrams. High 
redshift galaxies are often offset from the local star form¬ 
ing locus toward that of AGN (Mannucci et al. 2010; 
Kewley et al. 2013; Shapley et al. 2014; Steidel et al. 
2014). Spatially resolved data have shown that in some 
galaxies this is due to a combination of star formation 
and a weak active galactic nucleus (Wright et al. 2010; 
Trump et al. 2011), while others show pure star forma¬ 
tion with offsets seen in individual giant star forming re¬ 
gions (Jones et al. 2013b). Alternatively, recent work has 
shown that N-based indicators could be biased tracers of 
Oxygen (Amorm et al. 2010; Andrews & Martini 2013; 
Masters et al. 2014; Shapley et al. 2015). Distinguishing 
these options is critical for metallicity studies. 

Progress in measuring resolved emission line diagnos¬ 
tics at high redshift has been slow due to limitations 
in angular resolution and sensitivity. The most promis¬ 
ing efforts, targeting gravitationally lensed galaxies with 
adaptive optics (AO), have yielded only a few metallicity 
gradients to date (Jones et al. 2010, 2013a; Yuan et al. 
2013). Seeing-limited data have been obtained for larger 
samples (Queyrel et al. 2012; Troncoso et al. 2014), and 
they are in tension with AO-based results, highlighting 
the need for high spatial resolution. The sensitivity of 
field HST surveys (e.g. WISP, 3D-HST, PEARS) is usu¬ 
ally sufficient to study only the very brightest systems 
and their angular resolution does not benefit from lens¬ 
ing magnification. 

The mass-metallicity-SER relation has been measured 
by numerous groups in large samples extending to z 
3.5 for the most massive {M^ > 10^-^ Mq) and luminous 
sources (e.g., Mannucci et al. 2010). Below z ^ 1, faint 
sources can be probed with optical spectra (Henry et al. 
2013). At higher redshifts, lensed galaxy surveys have 
extended this work to lower and SER (Belli et al. 
2013), but are limited by small samples (i.e., only ^ 15 
high redshift galaxies with < 10^ Mq). The vast 
majority of these studies rely on integrated spectra and 
therefore cannot distinguish the cause of offsets in the 
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BPT diagram. 

GLASS addresses these issues by collecting an unprece¬ 
dented sample of emission line galaxies at 0.65 < z < 3.3, 
with spatially resolved information for a significant frac¬ 
tion. The questions raised above are addressed via the 
integrated flux and spatial distribution of multiple emis¬ 
sion line diagnostics which are sensitive to the gas metal- 
licity and ionizing source. In particular the strong lines 
[O II], [O III], and H/d are all observed for galaxies at 
1.3 < z < 2.3 providing good diagnostics of both metal- 
licity and nuclear activity (for example via the ’’blue” di¬ 
agnostic diagram; Lamareille 2010). [S II] and H(a+[N II] 
are also valuable diagnostics available at z < 1.5, while 
metallicities can still be estimated for the strongest line 
emitters at 2.3 < z < 3.3 via [O II] and [Ne III]. As 
shown by Jones et al. (2015), the quality of the GLASS 
data enables accurate measurements of metallicity gra¬ 
dients, and of the mass-metallicity-SFR relation down 
to M*=10^ Mq at z = 2. This progress is made pos¬ 
sible by the observing strategy and lensing magnifica¬ 
tion which deliver an order of magnitude improvement 
in flux sensitivity and of 3-4 in spatial resolution over 
previous HST grism surveys in blank fields (e.g., Henry 
et al. 2013). In general, the broad wavelength coverage of 
GLASS enables the comparison of metallicity estimates 
based on different features, so that one can test, for ex¬ 
ample, whether the [N II] based estimates are biased with 
respect to [S II]. 

2.3. Key Science Driver 3: Environmental Dependence 
on Galaxy Evolution 

The epoch 0 < z < 1 is one of rapid decline in the 
global star formation rate (Madau et al. 1996; Lilly et al. 
1996; Madau & Dickinson 2014), but clusters experience 
an evolution in star formation activity over this time that 
is even stronger than the general field. Identifying the 
processes that trigger and terminate star formation in 
cluster galaxies (Butcher & Oemler 1984; Dressier et al. 
1999; Poggianti et al. 1999; Dressier et al. 2013), and 
contrasting them to those operating in the field (Cooper 
et al. 2008; Vulcani et al. 2010; Muzzin et al. 2012; Oem¬ 
ler et al. 2013) is key to understanding the causes of 
the general decline. For what galaxy types and masses 
are environmental effects driving the quenching of star 
formation, and what is instead internally driven (Peng 
et al. 2010)? Is there one process that dominates over 
all environments or do some play a larger role in driv¬ 
ing galaxy evolution in the field than they do in dense 
environments? Massive clusters are the sites where the 
effects of environmental processes, such as ram pressure 
stripping, are expected to be strongest. Hence, studies 
of massive clusters are needed to answer these questions. 

Each of the processes that have been proposed to 
quench star formation in galaxies should leave a different 
signature on the spatial distribution of the star forma¬ 
tion activity within the galaxy. As an example, the loss 
of the galaxy hot gas halo would lead to a suppression of 
star formation over the whole disk, while ram-pressure 
stripping should leave a recognizable pattern of star for¬ 
mation with truncated Ha disks smaller than the undis¬ 
turbed stellar disk (see Yagi et al. 2015, for an example 
of ram pressure stripping). What is needed is thus an 
unbiased census of the star formation activity in a sta¬ 
tistically significant sample of clusters, yielding both the 


total star formation rates and the spatial distribution of 
star-forming regions within galaxies. A most interesting 
epoch for such a study is in the range z ^ 0.3 — 1, where 
most of the evolution takes place. 

Ha is considered the best optical indicator of the cur¬ 
rent star formation rate, being much less affected by dust 
extinction and metallicity effects than [O II] (Kennicutt 
1998; Moustakas et al. 2006). For this reason, a num¬ 
ber of Ha surveys up to z ^ 1 have been undertaken in 
the field using narrow-band imaging (e.g., Sobral et al. 
2013) and with WFG3 grism observations (e.g., Atek 
et al. 2010; Straughn et al. 2011). In clusters, narrow- 
band Ha studies are available for just a few systems at 
z = 0.3—1 (Kodama et al. 2004; Finn et al. 2005; Koyama 
et al. 2011) and a few other higher-z overdense regions 
(Koyama et al. 2013). These studies provide integrated 
Ha fluxes, and no spatial distribution information. The 
power of spatially resolved information is exemplified by 
a study of Virgo which suggested that in the local Uni¬ 
verse stripping of gas is the main mechanism for quench¬ 
ing star formation rates in clusters (Koopmann & Ken¬ 
ney 2004). 

GLASS is designed to measure Ha fluxes for all clus¬ 
ter star-forming galaxies in the central 0.6-0.9 Mpc of 10 
clusters at z = 0.31 — 0.69. Note that Ha and [N II] are 
blended at the resolution of G102; this effect can be taken 
into account when comparing with other samples, or cor¬ 
rected based on higher spectral resolution data (Trump 
et al. 2011, 2013; Fumagalli et al. 2012). GLASS’S target 
sensitivity is an order of magnitude improvement over 
previous studies of emission lines in cluster galaxies (So¬ 
bral et al. 2013). GLASS is designed to reach SFR limits 
comparable to the deepest narrow band studies, but for 
all objects and with spatially resolved information (see 
Vulcani et al. 2015, in preparation, for more details). 
Thus, GLASS is the first unbiased source of spatially 
resolved star formation within cluster galaxies at these 
redshifts. The combination of GLASS spectroscopy and 
HST photometry from GLASH/HFF is used to identify 
uniquely the line as Ha, as well as measure morphology, 
spatially resolved colors and stellar masses for virtually 
all objects detected in emission. 

Parallel AGS G800L Grism spectroscopy in two offset 
fields (at 2-3 Mpc, i.e. 1-3 virial radii; 2 position angles; 
see Figure 3) for each cluster, will cover the cluster infall 
regions, a key region for some environmental processes 
(e.g., Treu et al. 2003). In these regions H/? (and [O II] 
at z > 0.6) cab be used as star formation indicators, 
while [O III] further aids in redshift identification. By 
combining cluster cores, infall regions, and field, GLASS 
enables the study of the spatially resolved star forma¬ 
tion across a wide range of environments. Given the low 
number density of cluster members at such clustercentric 
distances, the AGS grism is the ideal method to build 
complete samples of star-forming galaxies. The parallel 
AGS data are not part of this first data release and will 
be discussed in detail in future publications. 

2.4. Additional Science Drivers 

GLASS spectroscopy and imaging is suitable for a wide 
range of applications. Although the survey design is 
driven by the three science cases described above, two 
additional science drivers were taken into account when 
possible. They are briefly summarized below. 
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a) Imaging 


b) GLASS Slitless spectroscopy 


Ebeling - ACS/WFC 
Siana - ACS/WFC 
CLASH - ACS/WFC 
CLASH - WFC3/IR 
HFF-ACS/WFC 
HFF - WFC3/IR 


GLASS - F105W/G102, F140W/G141 
GLASS - F814W/G800L 


Fig. 3.— Observational layout of the GLASS spectroscopy in the context of existing imaging data for MACSJ0717.5+3745. The WFC3 
grism spectroscopy is aimed at the cluster center overlapping with deep multiband imaging from CLASH and/or HFF, while the two parallel 
fields are approximately 90 degrees apart, with one of them coinciding with the CLASH/HFF parallels. 


2.4.1. Luminous and dark matter in clusters of galaxies 

Measuring the distribution of matter in clusters of 
galaxies is a powerful way to address a number of funda¬ 
mental astrophysical questions, ranging from cosmology 
to the interaction between supermassive black holes and 
their environment. What is the self-interaction cross sec¬ 
tion of dark matter? Is dark matter cold? What is the 
net effect of energy feedback from black hole accretion 
on dark matter halos? 

By modeling the gravitational lensing distortion of 
background galaxies one can obtain precise and accu¬ 
rate two-dimensional maps of the mass distribution in 
the clusters themselves. Considerable effort has gone and 
is going into developing modeling techniques that exploit 
all the available information, and deep HST imaging data 
have allowed researchers to identify more than one hun¬ 
dred multiple images of background sources per clusters. 
However the vast majority of these multiple images do 
not have spectroscopic redshift identification, and photo¬ 
metric redshift must be used to transform their positions 
into constraints on the gravitational potential of the clus¬ 
ters. Spectroscopic redshifts help enourmously by reduc¬ 
ing the effective uncertainty of each potential measure¬ 
ment to below 1% and by eliminating mis-identifications 
and catastrophic errors on photo-z which are common for 
extremely faint galaxies with unusual shapes. One of the 
ancillary drivers of GLASS is to measure spectroscopic 
redshifts for as many faint sources as possible in order to 
provide more extended and cleaner catalogs of multiple 
images to be used as input to lens models. First exam¬ 
ples of this application of GLASS are given by Schmidt 
et al. (2014b) and Wang et al. (2015). 

2.4.2. Supernovae 

Gravitational telescopes magnify all sources in the 
background, including transient ones. Thus, in princi¬ 


ple gravitational telescopes can be used to detect and 
take spectra of supernovae that are fainter and more dis¬ 
tant than otherwise possible (Brammer et al. 2014), thus 
extending the look back time over which supernovae and 
supernova cosmology can be studied (e.g. Sullivan et al. 
2000; Patel et al. 2014). If the lensed supernovae are 
of type la, knowledge of their absolute magnitude pro¬ 
vides an opportunity to break the mass-sheet degeneracy 
(Zitrin et al. 2014; Nordin et al. 2014) (Rodney et al. 
2015). 

In rare circumstances, if the supernovae happens to ex¬ 
plode within the cluster’s caustics, the supernovae can 
appear multiply imaged to the observer. The differ¬ 
ence in arrival time between the images is measurable 
and can be used as a cosmic ruler to measure distances 
and hence cosmology (Refsdal 1964; Holz 2001; Bolton & 
Buries 2003), provided the gravitational potential of the 
cluster can be sufficiently well constrained. This tech¬ 
nique has been demonstrated to be very powerful using 
quasars as the variable lensed source. Typically, the de¬ 
flectors are galaxies, for which the gravitational potential 
can be sufficiently well constrained by extended gravi¬ 
tational arcs and stellar kinematics (Treu & Koopmans 
2002; Suyu et al. 2014). Multiply imaged supernovae are 
significantly more difficult to find, and the first example 
has been discovered only recently in data taken by the 
GLASS survey (Kelly et al. 2015). 

This area of science has been done in close coordina¬ 
tion with the FrontierSN Program (GOs 13386 and GO 
13790; PI: Rodney) for the HFF clusters, ensuring uni¬ 
form acquisition and analysis of the data. 

3. SAMPLE SELECTION AND PROPERTIES 

The number and choice of clusters depend on a num¬ 
ber of factors, starting from the steepness of the number 
counts of the target population. The epoch of reioniza¬ 
tion science driver was used to select the optimal strat- 
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egy. Figure 2 shows the expected number of Lya emit¬ 
ters as a function of total number of orbits, using the 
model described in Section 2.1.1, plus the actual sensi¬ 
tivity measured for the GLASS grism data^. At z ^ 6—8, 
once a sufficient depth per cluster is reached, the num¬ 
ber of sources is maximized by observing more clusters 
rather than going deeper. Observing multiple clusters 
also minimizes cosmic variance (for z > 6 galaxies < 10% 
for 10 clusters, compared to 30% for a single WFC3 field 
Trent! & Stiavelli 2008). Observing only 4 clusters (e.g. 
the initial HFF targets) for a longer individual exposure 
reduces the number of detections because magnification 
becomes less effective. Of course, even though we chose 
the number of sources as metric, this is not the only rea¬ 
sonable choice. By going deeper on a smaller number of 
clusters one would have probed an intrinisically different 
population of galaxies. 

In order to find the 10 best clusters to target with 
GLASS we turned to the GLASH and HFF initiatives 
for the unique data that are being accumulated by a va¬ 
riety of groups with HST and other facilities. This vast 
array of data makes the identification of emission lines 
much easier than for clusters with just a few HST fil¬ 
ters, e.g. allowing for the use of the dropout technique 
to confirm that a single emission line is indeed Lya, and 
of photo-zs to resolve cases with multiple line identifica¬ 
tions. We further pruned the list of GLASS and GLASH 
clusters, using cluster redshift as an additional criterion. 
In fact, clusters at z ^ 0.3 —0.8 are the best gravitational 
telescopes for two reasons: i) the size of the critical lines 
is well matched to the WFG3 field of view; ii) cluster 
galaxies appear smaller and fainter than in lower red- 
shift clusters, thus minimizing the foreground contam¬ 
ination for background sources. The final requirement 
is given by the ability to observe Ha in cluster galaxies 
{z = 0.22 — 1.59). Based on these criteria we selected the 
sample of ten clusters at ^ = 0.31 — 0.69 listed in Table 1. 
The table lists the main (but not only) sources of HST 
imaging and Spitzer Imaging, as well as basic features 
like redshift. X-ray temperature and weak lensing virial 
mass. Much ancillary data is available for these clusters 
at virtually all wavelengths. 

4. OBSERVATIONS AND DATA REDUCTION 

The two grisms G102 and G141 provide continuous 
wavelength coverage in the range 0.81 — 1.69/im (defined 
as the range where the sensitivity is > 20% of the peak), 
thus enabling the detection of Lya above z = 5.7, metal- 
licity gradients over a range of redshifts, and Ha for the 
cluster members (see Figure 4). 

The depth of 10 orbits in G102 was driven by the need 
to reach a sensitivity sufficient to probe the Lya luminos¬ 
ity function at 2 : > 6. The depth of four orbits in G141 
was chosen to reach approximately uniform sensitivity 
across the two grisms, which is appropriate to measure 
typical line ratios required for metallicity determinations 
(Atek et al. 2010; Henry et al. 2013). In this section 
we detail the observational strategy (4.1) and the data 
reduction process (4.2). 

^ As it will be shown later in the paper, the sensitivity to Lyo; 
varies very rapidly between 2 : = 5.5 and 6.5. We are taking as a 
reference the value at 2 ; = 6, although most objects will be detected 
above 2 ; > 6 where the sensitivity is highest. 



Fig. 4.— Illustration of the wavelength coverage provided by the 
two WFC3 grism (in red and blue, respectively GI4I and GI02). 
The continuous wavelength coverage between ~ 0.81 — I.69/rm en¬ 
ables all three main science drivers. The transmission curves of 
the filters used by the CLASH imaging program are shown as light 
black curves for comparison. Four example spectra, one for each 
of the diagnostic features shown in the caption, are overplotted for 
illustration. 

4.1. Phase II design 

The main pointing was chosen to overlap with existing 
or planned HST imaging observations. The observations 
of each cluster were then divided into two sets of visits 
at different roll angles, to help dealing with the effects of 
contamination in the crowded cluster fields. One position 
angle was chosen so that its ACS parallel field would land 
in a HFF or CLASH parallel field, while the second roll 
angle was chosen to be at 90 ± 10 degrees from the first, 
so as to resolve contamination from overlapping spectra, 
while maximizing the overlap between the two WFC3 
grism pointings. 

A dither strategy analogous to that followed by the 3D- 
HST survey shown in Figure 3 of Brammer et al. (2012) 
has been adopted, in order to maximize spatial and spec¬ 
tral resolution and defect removal. For each visit at least 
4 subexposures were taken for each grism (typically half 
orbit each) with semi-integer pixel offsets, so that they 
could be combined by interlacing as described in the next 
section. 

A paired direct image exposure was taken with each 
grism exposure (F105W or F140W for G102 or G141, 
respectively) without offsetting the telescope for image 
alignment and spectral calibration. For visits that con¬ 
tained only G102 spectroscopy, some of the pre-imaging 
was obtained in F140W in order to ensure uniform filter 
coverage of transient events like Supernovae. This choice 
had virtually no impact of the depth of the F105W ex¬ 
posures or accuracy of the calibration of the G102 spec¬ 
troscopy while at the same time assisting for the ancillary 
science case. Furthermore, after aligning the grism vis¬ 
its to existing observations from the GLASH and HFF 
programs, the deeper images from those programs can 
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TABLE 1 

Sample properties 


cluster 

RA 

(J2000) 

DEC 

(J2000) 

z 

HST imaging 

Spitzer imaging 

Lx 

(lO^^erg s“^) 

M500 

First Release 

A2744 

00:14:23.4 

-30:23:26 

0.308 

HFFl 

SFF 

15.28+0.39 

17.6+2.3 

F15 

A370 

02:39:52.8 

-01:34:36 

0.375 

HFF3 

SFF 

8.56+0.37 

11.7+2.1 

W16 

M0416.1-2403 

04:16:09.4 

-24:04:04 

0.420 

CLASH/HFFl 

SURFSUP 

8.11+0.50 

9.1+2.0 

W16 

M0717.5+3745 

07:17:31.6 

+37:45:18 

0.548 

GLASH/HFF2 

SFF 

24.99+0.92 

24.9+2.7 

This paper 

M0744.9+3927 

07:44:52.8 

+39:27:24 

0.686 

GLASH 

SURFSUP 

18.94+0.61 

12.5+1.6 

W16 

M1149.6+2223 

11:49:35.9 

+22:23:55 

0.544 

GLASH/HFF2 

SURFSUP 

17.25+0.68 

18.7+3.0 

F15 

RXJ1347.5-1145 

13:47:30.6 

-11:45:10 

0.451 

GLASH 

SURFSUP 

47.33+1.2 

21.7+3.0 

S15 

M1423.8+2404 

14:23:47.8 

+24:04:40 

0.545 

GLASH 

SURFSUP 

13.96+0.52 

6.64+0.88 

S15 

M2129.4-0741 

21:29:26.1 

-07:41:29 

0.570 

GLASH 

SURFSUP 

13.69+0.57 

10.6+1.4 

S15 

RXJ2248 

22:48:44.3 

-44:31:36 

0.348 

GLASH/HFF3 

SFF 

30.81+1.57 

22.5+3.3 

S15 


Note. — For each cluster we list J2000 coordinates, redshift, the main sources of HST and SST imaging, as well as M 500 from 
Mgas and X-ray luminosity (from Mantz et al. 2010). The last column lists the target date for the first public data release (F=Fall, 
S=Summer, W=Winter). 


be used as the reference for the spectral extractions and 
modeling. 

The effective exposure times for MACSJ0717.5+3745 
in G102, G141, F105W and F140W at the two different 
position angles are listed in Table 2 and are typical of 
all GLASS clusters. The notional depth of the G102 and 
G141 grism exposures are 10 and 4 orbits respectively, 
including alignment images and overheads. 

The depth of the spectroscopic data is illustrated in 
Figure 5 which gives the sensitivity for an unresolved 
emission line in a clean part of the field of view. The 
sensitivity estimates were estimated for a sample of non¬ 
detections z '> 7 compact dropouts from the first 6 
GLASS clusters as presented by Schmidt et al. (2015). 
We used extraction apertures of 5 (spatial) by 3 (spec¬ 
tral) native pixels, which corresponds to apertures of 
^0.6 "x 100A similar to what was used by Schmidt et al. 
(2014b). For more extended emission and larger objects, 
a larger aperture of, e.g., 10 (spatial) by 6 (spectral) 
native pixels, might be more representative of the sensi¬ 
tivity estimates. Such an aperture will increase the noise 
level by a factor of 2. These noise levels are in fair agree¬ 
ment with previously published sensitivities of the HST 
NIR grisms (Brammer et al. 2012; Trump et al. 2014). 
We note that the depth of the spectroscopic data varies 
significantly from cluster to cluster owing to the variable 
background level, due to the combined zodial light and 
atmospheric emission. The depth of the imaging data 
is comparable to that of GLASH. All imaging data ob¬ 
tained as part of the GLASS program are combined with 
imaging data taken as part of the GLASH and HFF pro¬ 
grams. 

4.2. Data reduction 

As noted in Section 4.1, the GLASS observations are 
designed to comply with the 3D-HST observing strategy 
and were processed with an updated version of the 3D- 
HST reduction pipeline^ described by Brammer et al. 
(2012). The updated pipeline combines the individ¬ 
ual exposures into mosaics using AstroDrizzle (Gonzaga 
et al. 2012), replacing the MultiDrizzle package (Koeke- 
moer et al. 2003) used in earlier versions of the pipeline. 
The individual exposures and visits are aligned using 
tweakreg and grism sky backgrounds are subtracted us- 

^ http://code.google.com/p/threedhst/ 



Fig. 5.— The observed Icr noise level (not accounted for lens- 
ing magnification) for a sample of non-detections high- 2 ; dropouts 
from the first 6 GLASS clusters (see Schmidt et al. (2015) for de¬ 
tails). An extraction parameter of approximately O.O'+IOOA run¬ 
ning over the full G102+G141 wavelength range was used. Due to 
significantly varying contamination for individual objects between 
the two observed GLASS position angles, the noise level is shown 
for spectra extracted from a single P.A. Hence, the displayed limits 
correspond to the depth obtained from half of the GLASS data on 
each object. Gombining spectra at the two P.A. for contamination 
free wavelength ranges will lower the noise level by a factor of \/2- 
For more extended galaxies, more representative noise levels can 
be obtained using a larger extraction aperture, e.g., of 1.2"x200A, 
which will decrease the sensitivity by a factor 2. 

ing master sky images as described by Brammer (2015, in 
preparation), Kiimmel et al. (2011), and Brammer et al. 
(2012). The direct images were sky subtracted by fitting 
a 2nd order polynomial to each of the source-subtracted 
exposures. Each exposure is then interlaced to a final 
image with a pixel size of ~ 0706x 12(22)A for the 

G102(G141) grisms. Before sky-subtraction and interlac¬ 
ing each individual exposure was checked and corrected 
for elevated backgrounds due to the He Earth-glow de¬ 
scribed by Brammer et al. (2014). The pipeline is re-run 
in a way that does not flag bad reads as cosmic rays, as 
discussed by Brammer (2015, in preparation), using the 
script available here https : //github. com/gbrammer/ 
wfc3/blob/master/reprocess_wfc3 .py. The final in¬ 
terlaced and sky-subtracted mosaics of the G102 and 
G141 grism for the two GLASS position angles are shown 
in the 4 bottom panels of Eigure 6. 

Erom these final mosaics, the spectra of each individ- 
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Fig. 6 .— The MACSJ0717.5+3745 GLASS fields-of-view. The color composite image (top) is based on the CLASH imaging with blue, 
green and red channels as noted on the right. The four bottom panels show the final interlaced sky-subtracted G102 (left) and G141 (right) 
grism mosaics at the position angle 20 (top; green polygons) and 280 (bottom; magenta polygons) degrees. The individual spectra are 
extracted from these mosaics based on the extent of the corresponding object in the direct image mosaics (not shown here). 
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TABLE 2 

MACSJ0717.5+3745 Exposure Times 


Filter 

P.A.=20 

Uxp/[s] 

P.A.=280 

Uxp/[s] 

Total 

Uxp/[s] 

GI02“ 

9629 

10829 

20459 

GI4I“ 

3812 

4312 

8123 

FI05W 

1979 

1979 

3959 

FI40W 

712 

712 

1423 


Note. — The table summarizes the 
exposure times of the GLASS data for 
MACSJ0717.5+3745. These are typical 
of the GLASS dataset for all 10 clusters. 

The position angles listed correspond to the 
HST keyword PA_V3. The P.A. of the y- 
axis corresponds to PA_V3+44.69. Addi¬ 
tional imaging data are available from the 
HST archive. 

“Exposure times after correction for He 
Earth-glow (see Section 4.2) 

ual object are extracted by predicting the position and 
extent of each two-dimensional spectrum based on the 
SExtractor (Bertin & Arnouts 1996) segmentation of 
the corresponding direct images. As this is done for ev¬ 
ery single object, the contamination, i.e., the dispersed 
light from neighboring objects in the direct image field- 
of-view, can be estimated and accounted for. We note 
that a complete description of the 3D-HST image prepa¬ 
ration pipeline, spectral extractions, and spectral fitting, 
will be provided by Momcheva et al. (2015, in prep) 

5. DATA QUALITY, VISUAL INSPEGTION, AND 
GONSTRUGTION OF THE REDSHIFT GATALOG 

Given the richness and complexity of the data set, we 
find that visual inspection of a magnitude limited subset 
of the data is a very useful addition to the automatic 
processing. It was decided to carry out visual inspec¬ 
tion of the redshift determinations as well, in order to 
flag catastrophic failures in the automated redshift de¬ 
termination. The automated redshift determination and 
subsequent visual inspection are described in Section 5.2. 
In addition to the complete magnitude limited catalogs 
we also carried out a visual search for isolated emission 
lines through the entire dataset. This search yielded ad¬ 
ditional redshifts and is described in Section 5.3. 

5.1. Visual inspection strategy 

Visual inspection carried out with the publicly avail¬ 
able Graphic User Interface GIG (GLASS Inspection 
GUI) described in Appendix A. GIG provides a con¬ 
venient and efficient way to browse the dataset. The 
inspection is aimed at assessing the quality of the data 
and it is therefore independent of any additional analysis 
like photometric redshifts. 

The main goals of the GLASS team visual inspection 
were the following: 

1. Identify and flag catastrophic failures of the 
pipeline. Even though these are rare in numbers, 
they can create substantial problems in subsequent 
analysis if they are not properly flagged. Exam¬ 
ples include systematic errors in the reconstruction 
of the contamination spectra and artifacts due to 
edge effects. 


2. Assess the degree of contamination in the spectra. 
The degree of contamination of a spectrum depends 
on a combination of factors, for example the rela¬ 
tive surface brightness of the main object to that of 
the contaminants, their distance on the sky and on 
the presence of spectral features. Although some of 
these quantities can be calculated by the pipeline 
it is useful to provide an additional comprehensive 
human assessment of the degree of contamination. 
After significant experiments we decided to divide 
the spectra in three classes of contamination mild, 
moderate and severe. Roughly speaking, the three 
categories refer to spectra where the area of the 
detector with contamination sufficient to affect the 
target spectrum is considered to be respectively 
< 10%, 10 — 40% > 40% although the assessment 
takes into account all the factors listed above. As 
shown in Eigure 7, approximately 40% of all spec¬ 
tra have mild contamination. However, thanks to 
two-PA strategy, there is at least one clean spec¬ 
trum for more than 60% of all the objects in the 
catalog. 

3. Elag and identify strong emission lines and the 
presence of a continuum. During the inspection of 
the spectra a box is ticked to denote the detection 
of a continuum, and the wavelengths of identified 
lines are marked. This information can be used 
in the next step of determining redshifts, allowing 
one, for example, to visually inspect only the spec¬ 
tra where continuum or emission lines are present. 

4. Note additional features. The user can also provide 
additional comments, or flag the object as being a 
star or defect etc. 

In practice, in order to provide some safeguard against 
the inevitable subjectivity associated with visual inspec¬ 
tion, each cluster catalog is inspected by at least two 
team members. Eor MAGSJ0717.5+3745 co-authors 
T.T. and B.V. visually classified the entire dataset of 
1151 objects down to magnitude TLab < 24 (E140W). 
This deep classification was conducted on the first clus¬ 
ter for exploration purposes. Eor the remaining clusters 
full visual inspection will be limited to Hab = 23, plus 
a search for emission lines in the fainter objects). Eor 
reference the distribution of photometric redshifts, taken 
from the GLASH catalog, of the parent sample, and the 
sample limited to H < 23 is shown in Eigure 8. The two 
inspection catalogs were merged by averaging numeric 
flags, i.e. if one inspector assigned mild (0) contamina¬ 
tion and the other assigned severe (1) contamination, the 
entry in the combined catalog is 0.5. 

We caution that the inspection performed by the team 
is of general purpose, and even though it should provide 
a useful guidance for anyone interested in the GLASS 
data, it may not suffice for very specific applications. 
Eor example, a dedicated re-inspection of photometri¬ 
cally selected high redshift galaxies was carried out by 
the team while looking for faint Lya emission (Schmidt 
et al. 2015). 
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TABLE 3 
Redshift catalog 


ID* 

RA 

(J2000) 

DEC 

(J2000) 

z 

quality 

mult 

Notes 

036 

109.39293627 

37.77049774 

0.0000 

4.0 

0.0 

star 

055 

I09.398I0224 

37.770416861 

0.8785 

3.5 

0.0 


057 

I09.4003863I 

37.770271353 

0.0000 

3.0 

0.0 

star 

145 

I09.402300I6 

37.766627582 

1.3575 

2.0 

0.0 


173 

109.39848457 

37.766442326 

0.5560 

3.0 

0.5 


234 

109.39639585 

37.764684116 

0.5490 

4.0 

0.0 


236 

I09.393665I9 

37.764470231 

0.3900 

4.0 

0.0 


272 

109.39400607 

37.764400472 

0.3895 

4.0 

0.0 


273 

109.40514662 

37.764210184 

0.0000 

4.0 

0.0 

star 

299 

109.39108321 

37.763095453 

1.9000 

1.5 

0.0 


307 

109.39619713 

37.76335404 

2.5400 

4.0 

1.0 

image 13.3 


Note. — First entries of the redshift catalog. The full catalog is given in its 
entirety in the electronic edition. The column “quality” contains the quality 
flag as described in the text. The column “mult” is set to one when a single 
line is detected and there are multiple possible redshift interpretations. The 
column “note” lists special comments about the object, e.g. if the object is 
part of a known multiply image system. * In the flnal catalog ids are formatted 
as either Oxxxx or 9xxxx. The IDs with a leading 9 refers to a reduction, which 
used more aggressive SExtractor de-blending and detection parameters. This 
was done to accommodate the detection of objects near bright objects, which 
are not assigned individual IDs when using the CLASH SExtractor parameter 
file (IDs with a leading 0). 



Contamination level 

Fig. 7.— Contamination statistics for the MACSJ07I7.5+3745 
catalog. The distribution of contamination levels (0=mild, Immod¬ 
erate, 2msevere) assigned by one classifier is shown by object and 
by spectrum, for the GI02 (blue) and GI4I (red) grisms. Approx¬ 
imately 40% of the spectra suffer from mild contamination and are 
therefore “clean”. However, thanks to the two position angle ob¬ 
serving strategy more than 60% of the objects have at least one 
clean spectrum. 

5.2. Redshift determination 

Redshift determination is performed in two steps. In 
the first step, templates are fit to each of the four avail¬ 
able grism spectra independently (G102 and G141 at two 
PAs each) to determine a posterior distribution function 
for the redshift. If available, photometric redshift distri¬ 
butions can be used as input priors to the grism fits in 
order to reduce computational time. Quite frequently. 



12 3 4 


z 

Fig. 8. — Distribution of photometric redshifts for the all the 
sources in the GLASH catalog of MAGSJ07I7.5+3745. The cluster 
redshift {z = 0.548) is indicated by a vertical dashed line. 

the information content of the four exposures is rather 
uneven, for example because some of the spectra are 
strongly contaminated, or maybe because emission lines 
fall in only one of the two grisms (G102 or G141) and 
no continuum is present in the other. Rather than com¬ 
bining blindly the four pdfs the team decided to undergo 
a step of visual inspection using the dedicated publicly 
available GLASS inspection GUI for redshifts (GIGz; de¬ 
scribed in Appendix A). With the help of GIGz the user 
can flag which grism fits are reliable or alternatively en¬ 
ter a redshift by hand if the redshift is misidentified by 
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z 


Fig. 9. — Redshift distribution of emission line sources identified 
from the WFC3 GLASS spectra in the field of MACSJ0717.5+3745. 
Note the peak at the cluster redshift 2 : = 0.548 (vertical dashed 
line). Stars are not included in the histogram and tally. 

the automatic procedure. As described in the appendix 
GIGz is fully interactive. This is carried out for all galax¬ 
ies down to a magnitude limit of H< 23, as well for ad¬ 
ditional samples as described in the next two sections. 

Using GIGz we assigned a quality flag to the redshift 
according to the following scheme (4=secure; 3=prob- 
able; 2=possible; l=tentative, but likely an artifact; 
0=no-z). These quality criteria take into account the 
signal to noise ratio of the detection, the possibility that 
the line is a contaminant, and the identification of the 
feature with a specific emission line. For example, Q=4 
indicates multiple emission have been detected with high 
signal-to-noise ratio and thus there is no doubt of the 
redshift identification; Q=3 indicates that either a single 
strong emission line is robustly detected and the redshift 
identification is supported by the photometric redshift, 
or that more than one feature is marginally detected; 
Q=2 indicates a single line detection of marginal qual¬ 
ity; Q=1 indicates that there is someting but it is most 
likely an artifact. We also ticked boxes to indicate secure 
identification of some of the more common and stronger 
lines. In some cases there is ambiguity about the iden¬ 
tification of a single emission lines. Those instances are 
marked during the inspection process. This procedure 
is carried out independently by each inspector and then 
their outputs are combined. 

5.3. Search for emission lines in faint sources 

In addition to the systematic visual inspection of the 
magnitude limited catalog, a search for emission lines was 
carried out by three of the authors (T.T., K.B.S., B.V.), 
based on visual inspection of all the jpegs of the two di¬ 
mensional spectra of objects with continuum fainter than 
the H < 23 limit. For each object, the authors inspected 
the spectra, the contamination model, and the contami¬ 


nation subtracted model. The two PAs and stacked spec¬ 
tra were simultaneously inspected together with postage 
stamp images of the system to allow for quick identi¬ 
fication of spurious contaminants, like zero order spec¬ 
tra from other sources. Gandidate emission lines were 
flagged for further inspection using GiGz, as described 
above. With the aim of providing a robust redshift cat¬ 
alog, as opposed to a complete redshift catalog, we took 
the conservative approach of considering as real lines only 
those identified by both human classifiers. This subjec¬ 
tive and quick procedure is sufficiently fast that allows a 
single investigator to look at all the over 20,000 spectra 
obtained by the survey and at the same time provides 
a first catalog of faint emission line objects. A more 
systematic search for emission lines using machine based 
methods (see, e.g., Maseda et al. 2015, in preparation) 
is left for future work. 

5.4. Search for multiple image redshifts 

Given the importance of measuring spectroscopic red- 
shifts for multiply imaged systems, all multiple image 
candidates compiled in the recent paper by (Diego et al. 
2014) were subject to an additional round of visual in¬ 
spection by two of the authors (T.T. and X.W.). We con¬ 
firm the redshifts of several previously known systems, 
including from our own GLASS survey (paper 0) and 
those measured by other teams (Ma et al. 2008; Limousin 
et al. 2012; Ebeling et al. 2014), and measure three new 
redshifts. One of them (image 5.1, GLASS ID=1108, 
z = 0.928) is classified as probable and differs substan¬ 
tially from the value z = 4.3 assumed by Diego et al. 
(2014). Two of the redshifts are tentative and should not 
be used until further confirmation (29.2, GLASS ID=378 
z = 1.73; 55.2, GLASS ID=82, z = 0.47, which would be 
in the foreground if confirmed). 

5.5. Properties of the spectroscopic sample and example 

spectra 

A histogram of redshifts based on emission lines de¬ 
tected in the WFG3 grism spectra is shown in Figure 9. 
Besides a clear peak associated with the cluster redshift, 
the distribution covers a broad range in redshifts. Some 
examples of spectra are shown in Figures 10 to Fig. 15. 
We show one example each of quality 4,3,2 redshifts, in¬ 
cluding multiply image candidates when possible. The 
examples are meant to illustrate the diversity of the data, 
including incomplete datasets, arising from edge effects 
or strong contamination. 

6 . PUBLICLY AVAILABLE DATA AND TOOLS 
6.1. Release plan 

For each cluster we plan three releases: 

1. Release 1. WFG3 images and spectra. In order 
to provide the most useful match to the photome¬ 
try provided by the GLASH collaboration for the 
first GLASS, the spectra are extracted based on 
photometric catalogs generated by running SEX- 
tractor (Bertin & Arnouts 1996) on the GLASS 
IR detection image using the same parameters as 
those adopted by GLASH (Goe, D. 2015, private 
communication; available at the GLASS website). 
All the spectra are visually inspected to magnitude 
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Fig. 10.— Spectra for GLASS ID ^1131 (Also known as arc 4.1). In each sub-figure, the two panels on top show the 1-dimensional 
spectra, where the observed fiux and contamination model are denoted by blue solid and red dashed lines respectively. The cyan shaded 
region represents the noise level. For the four panels directly underneath, the middle two display the interlaced 2-dimensional spectra 
whereas the bottom two have contamination subtracted. In the 1- and 2-dimensional spectra, the identified emission lines are denoted by 
vertical dashed lines in magenta and arrows in red respectively. The two panels on the left show the 2-dimensional postage stamp created 
from the HFF co-adds through drizzling (top) and the 1-dimensional collapsed image (bottom). Note that these two panels share the same 
x-axis along the grism dispersion direction. Some ancillary information is also shown in the upper left corner in each sub-figure, including 
the redshift quality fiag Q. 
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Fig. 11.— As Figure 10 for cluster member GLASS ID 431. 
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Fig. 12 .— As Figure 10 for GLASS ID 1402 (arc 12.1). The data at PA=280 suffer from strong contamination from the nearby galaxy. 
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Fig. 13.— As Figure 10 for cluster member GLASS ID 2334. 
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Fig. 14.— As Figure 10 for GLASS ID 90378 (arc 29.2). Data at PA=20 suffer from strong contamination from a nearby star. The data 
at PA=280 show low signal-to-noise ratio emission lines, thus justifying quality flag Q=2. 
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Fig. 15.— As Figure 10 for cluster member GLASS ID 1747. A low signal to noise ratio emission line is detected at both PA, and 
tentatively identified as Balmer Ha at the cluster reddshift. 
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F140W<23AB. The structure of the data and nam¬ 
ing of the files is identical to that produced for 3D- 
HST, with the addition of stack files that combine 
the two PAs into a single spectrum. More details 
are available from the 3D-HST paper and at the 
HST archive. In addition, a catalog is provided in 
electronic form. The first few entries are shown in 
Table 3 for guidance. 

2. Release 2. ACS images and spectra in the par¬ 
allel fields. The spectra are extracted based on 
catalogs generated from the GLASS ACS images 
themselves. Visual inspection will be performed to 
a magnitude limit to be defined. The data and cat¬ 
alog structure will be as similar as possible to that 
of the WFC3 data and will be described in a future 
publication. 

3. Release 3. WFC3 and ACS spectra. The spectra 
are extracted based on catalogs generated from the 
final images of the clusters, after the completion of 
the HFF campaigns. Spectra of newly identified 
objects are visually inspected and redshifts are re¬ 
calculated based on updated photometric redshifts. 

Barring the unexpected, the first release is scheduled 
to happen within approximately one year from the ac¬ 
quisition of the data (see Table 1; release will start 
after this paper is accepted). The second release is 
tentatively scheduled to follow by approximately 6-12 
months. The third release will depend on availability 
of other datasets and availability of funds. The pub¬ 
licly released data will be available through the Mikulski 
Archive for Space Telescopes (MAST) at URL https: 
//archive.stsci.edu/prepds/glass/. 

7. SUMMARY 

This paper provides an overview of the GLASS survey. 
After reviewing the scientific motivations that drove the 
survey design, the paper describes the data and data 
processing steps. We use MACSJ0717.5+3745, the first 
cluster observed by GLASS (and part of the HFF cam¬ 
paign), to illustrate the data set. The main results of 
this paper can be summarized as follows: 

1. The line flux sensitivity of GLASS reaches ^ 4 — 
5 • 10“^^ erg s“^ cm“^ (1"<^) foi* clean spectra of 
compact sources in the most sensitive wavelength 
range, for each of the two position angles. The ob¬ 
serving strategy, consisting of two orthogonal po¬ 
sition angles, is effective in reducing the effects of 
foreground contamination. Whereas ^40% of the 


individual spectra are classified as suffering from 
mild contamination (i.e. being virtually clean), 
more than 60% of the objects are observed as clean 
at at least one position angle. For those objects for 
which clean spectra are available at both position 
angles, the sensitivity reaches ~ 3 • 10“^^ erg s“^ 
cm“^ (l-<^)- 

2. The first data release for cluster 
MAGSJ0717.5+3745 is complete, and high 
level data products made publicly available. 
Spectra for 1151 galaxies down to magnitude 
Hab < 24 (F140W) have been visually inspected 
by members of our team to ensure quality control. 

3. A visual search for emission lines has been carried 
out through the entire dataset, including galaxies 
fainter than the inspection limit. In total, we mea¬ 
sure emission line redshifts for 139 extragalactic 
sources. The redshift catalog is made public in 
electronic format as part of the first data release. 

4. A dedicated search for redshift of candidate multi¬ 
ple imaged sources reveals three new redshifts, as 
well as confirming several previously known sys¬ 
tems. 

5. In addition to the high-level data products we 
make available two Graphic User Interfaces (GiG 
and GiGz) which allow for efficient browsing of 
the dataset, and interactive redshift determination. 
The GUIs are described in the appendix to this pa¬ 
per. 
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APPENDIX 


A. THE GLASS INSPECTION GUI 

In this appendix we describe the GLASS inspection GUI (GiG) and the GLASS inspection GUI for redshifts (GiGz) 
vl.O, released as part of the first public data release from GLASS described in this paper. This Python-based software 
provides a convenient and efficient way of inspecting and browsing the GLASS data products and is made publicly 
available for download at https://github.com/kasperschmidt/GLASSinspectionGUIs. A more detailed and con¬ 
tinuously updated description of the software package can be found in the GiG README which is also available at 
https://github.com/kasperschmidt/GLASSinspectionGUIs. 
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Fig. 16. — Overview of the GiG interface. The green text comments on the individual parts of the GUI. The software is publicly available 
for download at https://github.com/kasperschmidt/GLASSinspectionGUIs 


A.l. The GLASS inspection GUI (GiG) 

GiG is designed to inspect and browse the main products from the data reduction pipeline described in Section 4.2. 
Overview and descriptions of the products are available from the 3D-HST papers and the HST MAST archive. In brief, 
they include the extracted two-dimensional spectra, stacks of the spectra at the two separate position angles in 2D, 
extracted one-dimensional spectra, and various diagnostic plots. GiG (and GiGz described below) are self-contained 
in the Python script visualinspection.py and only depend on standard publicly available python packages (see the 
GiG README for details). A general overview of the interface of GiG v2.2 is shown in Eigure 16. GiG allows the 
user to quickly visualize all data products, including direct inspection of twodimensional fits files via interfacing with 
ds9. In addition to visualization, GiG allows the user to rate the contamination of the individual spectra at the two 
GLASS PAs, indicate the presence of a continuum, mark defects in the spectra or contamination models and mark 
emission lines including noting the redshift at which the lines were identified. 

The contamination levels are rated as mild, moderate, or severe corresponding to roughly <10%, 10-40% and >40 
flux contamination of the central region where the spectral trace of the object displayed is expected to be. As this 
rating is somewhat subjective an automatically generated estimated of the level of flux in the pixels ’belong’ to the 
spectrum is stored in the GiG output file. This contamination level is defined as the fraction of pixels in the fits- 
extension of the extracted 2D spectra containing the contamination model (masked by the object model extension) 
with > 10“^e/s/pixel. 

Apart from the manually set keywords and the automatic contamination estimate, the spectral coverage of each 
spectrum is also estimated and stored automatically in the output file. 
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Fig. 17.— Overview of the GiGz interface. The green text comments on the individual parts of the GUI. The software is publicly 
available for download at https://github.com/kasperschmidt/GLASSinspectionGUIs 


A.2. The GLASS inspection GUI for redshifts (GiGz) 

GiGz is developed for general inspection of the extracted ID spectra, the inspection of the redshift fits generated as 
described by (Brammer et al. 2012), and for manual redshift fitting of any detected emission lines. The default GUI 
window shown for GiGz vl.O in Figure 17 enables quality assessment of the redshift fits and flagging of particular lines 
identified for easy identification of, e.g., [O III] emitters in the data set. GiGz also includes an interactive plotting 
interface which can be controlled from the main window. This enables a convenient way of plotting and inspecting the 
extracted ID spectra, the redshift models, and any emission lines identified in the GiG inspection. Twodimensional 
fits files of the spectra can also be inspected through GiG via a ds9 interface. 
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